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UNIVERSAL UPPER AND LOWER BOUNDS ON ENERGY OF 

SPHERICAL DESIGNS 

P. G. BOYVALENKOV P. D. DRAGNEV D. P. HARDIN* *, E. B. SAFE*, 

AND M. M. STOYANOVA** 


Abstract. Linear programming (polynomial) techniques are used to obtain lower and 
upper bounds for the potential energy of spherical designs. This approach gives unified 
bounds that are valid for a large class of potential functions. Our lower bounds are 
optimal for absolutely monotone potentials in the sense that for the linear programming 
technique they cannot be improved by using polynomials of the same or lower degree. 
When additional information about the structure (upper and lower bounds for the 
inner products) of the designs is known, improvements on the bounds are obtained. 
Furthermore, we provide ‘test functions’ for determining when the linear programming 
lower bounds for energy can be improved utilizing higher degree polynomials. We also 
provide some asymptotic results for these energy bounds. 


1. Introduction 

Let be the unit sphere in M”. We refer to a finite set C C as a spherical 
code and, for a given (extended real-valued) function h : [—1,1] —)• [0,-t-oo], we consider 
the h-energy (or the potential energy) of C defined by 

(1) E{n,C]h):= ^ h{{x,y)), 

x,y£C,x^y 

where (x, y) denotes the inner product of x and y. At times we shall require h to be 
absolutely monotone or strictly absolutely monotone on [—1,1); i.e., its k-th derivative 
satisfies > 0 {h^^\t) > 0, resp.) for all fc > 0 and t G [—1,1). 
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A spherical r-design C C whose cardinality we denote by |C|, is a spherical code 

such that 

is the surface area measure) holds for all polynomials f{x) = f{xi,X 2 ,---,Xn) 
of total degree at most r. The maximal number r = t{C) such that (7 is a spherical 
r-design is called the strength of C. 

A commonly arising problem is to estimate the potential energy of certain sets of 
codes C. In this paper we address this problem for the class of spherical designs of fixed 
dimension, strength and cardinality. Denote by 

(2) £(n, N, r; h) := inf{£'(n, C; h) : IC] = A^, C d is a r-design} 

and 

(3) Z^(n, A", r;/i) := sup{£^(n, C;/i) : IC*! = A^, C C is a r-designj 

the minimum and the maximum possible /i-energy of a spherical r-design of N points 
on respectively. In this paper we derive lower bounds on C{n, N,t; h) and upper 

bounds on U{n,N,T-,h), which then define a strip where the energies of all designs of 
hxed dimension, strength and cardinality lie (see Theorem 13. 7 jl . 

Concerning lower bounds for energy, a general linear programming technique originally 
introduced by Delsarte m for investigating codes over finite fields (see also m) has been 
utilized by Yudin m, Kolushov and Yudin [30] and Andreev [2] (see also [D El ESI) to 
prove the optimality of certain spherical codes among all possible codes. In 2007, Cohn 
and Kumar augmented this technique by introducing the notion of conductivity to prove 
the universal optimality of shar^ codes. By universal optimality of a code we mean that 
among all codes of the same cardinality, it minimizes the energy for all absolutely mono¬ 
tone potentials. Essential to their method are certain quadrature rules associated with 
those sharp codes. Here, by combining the Delsarte method with quadrature formulas 
developed by Levenshtein, we derive lower bounds for general designs that are optimal 
for the linear programming technique when restricted to the use of polynomials of fixed 
maximal degree. 

Since the collection of spherical designs is a special subclass of codes for a given dimen¬ 
sion and cardinality, we expect that both larger lower bounds and lower upper bounds are 
possible when compared with minimal energy codes. Indeed we provide such examples 
in Section 5, where we discuss Mimura 2-designs (see Example 5.4). Our general upper 
bounds for A(n, A, r; h) are given in Theorems 13.61 and 13.71 in Section 3. 

We remark that upper bounds for the Riesz s-energy of well-separated spherical de¬ 
signs of asymptotically optimal cardinality on (that is, \C\ = 0{t‘^) as r —)• oo) were 


code is sharp if for some positive integer m, the code is a 2m — 1 design with at most m different 
values of the distance between distinct points in the code. 
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obtained by Hesse and Leopardi in [251 [26] . The existence of such well-separated de¬ 
signs of asymptotically optimal cardinality on (and more generally on has been 

established recently by Bondarenko, Radchenko, and Vyazovska in [8|. 

An outline of our paper is as follows. Some preliminaries are explained in Section 2, 
where we refer to results and techniques developed by Delsarte, Goethals and Seidel [21] 
and Levenshtein inaEniEn] that will be needed for the statements of our main results. 
The relationship between the Delsarte-Goethals-Seidel bounds for the minimum possible 
size of spherical designs of prescribed dimension and strength and the Levenshtein bounds 
for the maximum size of spherical codes of prescribed dimension and minimum distance 
will play a very important role in our investigation. Some results on the structure of 
designs of fixed dimension, strength and cardinality from mm are also discussed. 

In Section 3 we formulate two general results that provide the framework for obtaining 
lower and upper bounds for the energy of spherical designs. Theorem 13.II (lower bounds) 
is a slight modification of a known result (cf. [TT] Proposition 4.1], [9] Chapter 5]), but 
Theorem 13.61 (upper bounds) is new. 

Theorem [331 gives lower bounds that are optimal in the sense described in Theorem l3.5l 
- they cannot be improved by utilizing polynomials of the same or lower degree that 
satisfy the conditions of Theorem 13.11 Following Levenshtein [ 33 ] we call these bounds 
universal. 

Some of the lower bounds from Theorem [331 can be further improved by either restrict¬ 
ing the interval containing the inner products of even-strength designs (Theorems [331 and 
14.2p or by allowing polynomials of higher degree for odd-strength designs (Theorems 14.31 
and 14.41 and Corollary 14.5p . For the latter case, Theorem 14.31 provides necessary and 
sufficient conditions for the global optimality of the bounds from Theorem 13.41 

Section 5 is devoted to improving upper bounds for spherical designs utilizing restric¬ 
tions on their inner products. Some asymptotic results and numerical examples that 
illustrate these upper bounds are also included. 

Finally, in Section 6 we derive an asymptotic lower bound for the energy of spherical 
designs of fixed strength as the dimension and cardinality grow to infinity in certain 
relation. An example (Euclidean realization of the Kerdock codes [28]) illustrating the 
tightness of these asymptotic bounds is presented. 

2. Preliminaries 

The results to be presented in Sections 3-7 utilize the notations and fundamental facts 
described in the subsections below. 

2.1. Gegenbauer polynomials and the linear programming framework. For fixed 
dimension n, the Gegenbauer polynomials [ 30 ] are defined by = 1, = t 

and the three-term recurrence relation 

(i -|- n — 2) P^^\{t) = {2i + n — 2)tP^'^\t) — iP/r|(t) for i > 1. 
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We note that are orthogonal in [—1,1] with a weight (1 — and satisfy 

= 1 for all i and n. If f{t) G M[t] is a real polynomial of degree r, then f(t) can 
be uniquely expanded in terms of the Gegenbauer polynomials as 

r 

( 4 ) fit) = Y,hPt\t)- 

i=0 


We use the identity (see, for example, 1211 Corollary 3. 8], m Equation (1.7)], [Ml 
Equation (1.20)]) 

(5) ici/(i)+ ^ /((^,y)) = |ci2/o+E^f;(E^b(^) 

x,y£C,x^y 2=1 ^ j = l \x£C 


as a source of estimations by polynomial techniques. Here C C is a spherical code, 
/ is as in (|1]), {Yij{x) : j = 1, 2,... , r*} is an orthonormal basis of the space Harm(i) of 
homogeneous harmonic polynomials of degree i and r* = dim Harm(i). 

The Delsarte-Goethals-Seidel bound and the Levenshtein bound described in the next 
subsections are obtained after the sums on both sides of ([5]) are neglected for suitable 
polynomials. 


2.2. Delsarte-Goethals-Seidel bound for spherical designs. Denote by D{n, r) the 
Delsarte-Goethals-Seidel [2T] bound for spherical designs 


( 6 ) 


Bin, t) > D[n, r) 




if r = 2A; — 1, 


fn + k — kn + k — 2\ 

\ n — 1 / ~*~V n — 1 J ’ 


if T 


2k, 


where H(n,r) := minljCj : G C ^ is a spherical r-design}. This bound plays an 
important role in our (initially heuristic) choice of applications of Theorems 13.II and 13.61 
We shall utilize the values of the function D(n, r) to decide the degrees of the poly¬ 
nomials to be used for the bounding of C{n, N,t-, h) and Uin,N,T-,h). The rule is the 
following - if we have dimension n, strength r and cardinality N G {Din, r), D(n, r -|- 1)], 
then we use polynomials of degree r for the lower bounds and r or r — 1 (depending on 
the parity of r) for the upper bounds. 


2.3. Levenshtein bounds for spherical codes. We now formulate and discuss the 
Levenshtein bounds [Ml ESI IMl IM] on 

(7) ^in, s) := max{jCj : C C {x, y) < s for all x,y £ C,x ^ y}, 

the maximal possible cardinality of a spherical code on of prescribed maximal inner 
product s. 
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For every positive integer m we consider the intervals 

if m = 2A; — 1, 


■— * 


1,1 1,0 


. 1,0 . 1,1 

‘'k ’ ‘'k 


if m = 2k. 


Here = —1 by definition, a,b £ {0,1}, i > 1, is the largest zero of the Jacobi 

polynomial ^ ^ (t). The intervals {Xm}m=i well defined (see [Ml Lemmas 

5.29 and 5.30]) and therefore constitute a partition of X = [—1,1) into countably many 
closed subintervals with nonoverlapping interiors. 

For every s G X^, Levenshtein introduces a certain polynomial of degree 

m that satisfies all the conditions of the corresponding linear programming bounds for 
spherical codes and yields the bound 


(8) A{n, s) < < 


T _ (k+n-3\ r 2fc+n-3 _ (^) 1 r i, a T 

^ 2 k-i{n, s) J forsGX 2 fe-i, 

r . fk+n—2\ [2k+n—l (l+®)(Tfc ^(■*)“Tt+l(^)) 1 „ r- T 

I ) |,„J f” 0 6 


The function 
(9) 


X(re, s) := 


L2fc-i(n,s) forsGX2fc-i, 

L 2 k{n,s) for s G l 2 k 

is continuous in s. 

The connections between the Delsarte-Gothals-Seidel bounds and the Levenshtein 
bounds are given by the equalities 


( 10 ) 


L2fc-2(ra,4’^i) = X2fc-i(n,4’^^) = D{n,2k - 1) = 2^^^ ^ , 


(11) X2fc-i(n,4’°) = L 2 k{np]f) = D{n,2k) = 

occurring at the ends of the intervals Im- 


n + k — 1 
n — 1 


n + k — 2 
n — 1 


2.4. A useful quadrature. It follows from |34l Section 4] (see also [33l [T3]) that for 
every hxed (cardinality) N > D{n, 2k — 1) there exist uniquely determined real numbers 
— 1 < uq < cq < • • • < CKfc-i < 1 and po,pi,..., Pk-i, Pi > 0 for z = 0,1,..., A; — 1, such 
that the quadrature formula 


( 12 ) 


/o — 


r(n-l) 

2«-2f(^) 2 y_i 




k-l 


N 


i=0 


holds for every real polynomial f{t) of degree at most 2/c — 1. 
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The number a^-i = ak-i{N) in (1121) is the solution of the equation L{n,s) = N for 
s > Once it is determined, it is known (and easily verified) that the remaining 

Oj = ai{N), i = 0,1,... ,k — 2, are roots of the equation 

Pk{t)Pk-i{ak-i) - Pk{ak-i)Pk-i{t) = 0 , 

where Pi{t) = P^^ jg g, Jacobi polynomial of degree i. In fact, ai, i = 

0,1,... ,k — 1, are the roots of the Levenshtein polynomial |34l Eqs. (5.81) and (5.82)] 
(see also m Theorem 5.39]) 

f2k-t^\t )=(^ - «fc-i) 

i=0 

used for obtaining the bound L 2 k-i{n,s), s = ak-i, in (| 8 ]l. 

Similarly, for every fixed N > D{n,2k) there exist uniquely determined real numbers 
— 1 = /3o < /?! < • • • < < 1 and 70 , 71 ,..., 7 ^, 7 i > 0 for z = 0 , 1 ,..., fc, such that the 

quadrature formula 

k 

(13) = ^ + 

1=0 

is true for every real polynomial f{t) of degree at most 2k. The numbers /3j = j5i{N), 
i = 0,1,... ,k, are the roots of the Levenshtein polynomial 

fc-i 

Ak’^'^Ht) = it- Po)it - Pk) n(^ “ 

i=l 

(used for L 2 kin,s), s = fik, in dS])). 

As mentioned in the Introduction we always take into consideration where the cardi¬ 
nality N is located with respect to the Delsarte-Goethals-Seidel bound. It follows from 
the properties of the bounds D{n,T) and Lr{n,s) (see (fTUI) and (fTTI) ') that 

N G [D{n,T), D{n,T + 1)] <)=)> sGIt, 

where n, s and N are connected by the equality 

N = Lr{n, s). 

Therefore we can always associate N with the corresponding numbers: 

(14) ao,«!,..., ak-i,po, pi,-.-, Pk-i when N G (Zl(n, 2k - l),D{n, 2k)] 
or with 

(15) /3o,/3i,- ■ ■ ,Pk,lo,li,-■ ■,1k when N G {D{n,2k), D{n,2k -h 1)] . 
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2.5. Bounds on smallest and largest inner products of spherical designs. De¬ 
note 

(16) u{n,N,T) := sup{n(C') : C C is a r-design, \C\ = A^}, 
where u{C) := max{(x,y) : x,y G C,x y}, and 

(17) £{n,N,T) := inf{.^(C') : C C is a r-design, \C\ = iV}, 

where £{C) := min{(x, y) : x,y G C,x y}. 

For every n, r and cardinality N G [D(n, r), Il(n,r -|- 1)] non-trivial upper bounds 
on u{n,N,T) can be obtained (cf. |13( 111]). Similarly, for even t = 2k and cardinality 
N G [D{n, 2k),D(n, 2k + 1)) non-trivial lower bounds on £{n, N, 2k) are possible [131111] . 
We describe here explicitly the cases r = 2 and r = 4. 

In [131H] the quantities u{n, N, r) and i{n, N, r) are estimated by using the following 
equivalent definition of spherical designs, which is a consequence of (l5|): a spherical r- 
design C C is a spherical code such that for any point x G and any real 

polynomial of the form f{t) = YlJ=o the equality 

( 18 ) J]/((x,y)) = /o|C| 

xec 

holds. 


Lemma 2.1 ([H]). Let n > 3. 

(a) For every N G [D{n, 2), D(n, 3)] = [n -|- 1,2n] we have 


(19) 


N -2 

u{n, N, 2) <-1. 


n 


(b) For erery e [D(n, 4), D(n, 5)] = [n(ra-|-3)/2, n(n-|-1)] we have 


( 20 ) 


_ < 2(3+v'(~-l)l(n + 2)iV-3(„ + 3)l) _ 

n[n + 2) 


Proof. To prove (a) one can use the polynomial f{t) = yt + y - j in (HID with 

y being the midpoint of the geodesic arc between two of the closest points in C as in 
Theorem 3.2 from m- The bound in (b) is proved in the example following Theorem 
3.2 from [11]. □ 


Lemma 2.2 m)- Let n > 3. 

(a) For every N G [D{n, 2), D{n, 3)) 


[ra -|- 1, 2n) we have 
N 
n 


( 21 ) 


£{n,N,2) > 1 
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(b) For every N G [D{n,A),D{n,5)) 


( 22 ) 


i(n, N,A)>1-- 
n 


= [n(n + 3)/2, n{n + 1)) we have 


Proof. The bound in (a) can be proved by utilizing the polynomial f{t) = in (IlSp with 
y being the midpoint of the geodesic arc between points —x and z G C, where x, z is a 
pair of points of C with smallest inner product (see Theorem 3.3 from [H])- The bound 
in (b) is proved as in the example after Theorem 3.3 from |llj . □ 


We remark that different bounds on u{n, N, r) and i{n, N, r) can be obtained by a 
technique from m Sections 2 and 3]. For r > 4 such bounds are better in higher 
dimensions than those from Lemmas IQ and [221 More generally, when t = 2k we 
establish in Lemma 12.31 lower bounds on l{n, N, r) that will be used in Section 4 to 
establish Theorem \n\ 

Lemma 2.3 ([l3l Lemma 4.1]). Let N G {D{n, 2k), D{n, 2k + 1)) and 

/(t) = (t-/3i)2...(t-/3fe)2, 

where f3o,..., f3k are as in (USD. Let f and rj denote the smallest and largest roots, 
respectively, of f{t) = 7 oA^/(—1). Then ^ < i{n,N,2k) < u{n,N,2k) < p. 

Proof. For the convenience of the reader, we provide a proof of Lemma 12.31 Let (7 be a 
spherical design of strength r = 2/c on with jCj = N. Let x and y be distinct points 
in C. Since f{t) >0 for all t and / vanishes at j3i,... ,j3k, it follows from (fT 8 l) and (fT3]l 
that 

/(I) + y)) < ^ /((x, z)) = Nfo = /(I) + A^7o/(-1), 

zee 

and so f{{x,y)) < A^ 7 o/(—1). Since f{t) is strictly decreasing for t < f3i and strictly 
increasing for t > Pk, we must have < {x,y) < r/. □ 

Note that will produce a non-trivial bound ^ > — 1 if and only if 70 IV < 1 and the 
next assertion, implicit in |151 Section 4] (see also Remark 5.58 in [3l]), shows that this 
is indeed true for N G {D{n, 2k), D{n, 2k + 1)). 

Lemma 2.4. If N £ {D{n,2k), D{n,2k + 1)), then jqN G (0,1). Hence, i{n,N,2k) > 

- 1 . 


Proof. We have the formulas (Equation (5.113) in |34] 1 

Tk{s, 1 ) 


70 = 


7a;(- 1, -l)Tfc(s, 1) - Tk{-1, l)Tk{s, -1) 
and (the equation in the last line of page 488 in |34j l 

Tk{l, l)Tk{s, -1) - Tk{l,-l)Tk{s, 1) 


N = L2k{n,s) = 


Tk{s,-l) 
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A little algebra then shows that 


7o^ = 


T-A{s) 


r-i/A(s)’ 

where A(s) = Tk{s,l)/Ti^{s, —1) as in [15] and T = Tk{l,l)/Tk{l, —1). Moreover, we 


have 


A{s)=T 




from [m Lemma 5.24], where > 0 and < 0 for every s G (see 

Lemmas 5.29 and 5.30 in [34]). Therefore the signs of A(s) and T are opposite. We 
conclude that 

^ m + iA»)i ^ 

|r| + i/|A(»)| 

The ratio is decreasing in s in the interval (see [Ml Lemma 5.31]), 

i.e. lA(s)j is increasing in s G (^^1:’• Since = 0 and 1 for s = and 
respectively, we obtain that 'yoN increases from 0 to 1 when s increases from s = to 

ty. □ 


3. General lower and upper bounds 

The general framework of the linear programming bounds for the quantities £(n, N, r; h) 
and U{n,N,T;h) is given by the next two theorems. Theorem 13.11 is an adaptation of 
known results (cf. [41 1 117j ) to spherical designs, we are not aware of any prior use of 
linear programming techniques for obtaining upper bounds on the potential energy as in 
Theorem 13.61 

Theorem 3.1. Let n, N, r be positive integers with N > D{n,T) and let h : [—1,1] 

[0, +oo]. Suppose I is a subset of [—1,1) and f{t) = fiPi^^(t) ® polynomial 

sueh that 

(Al) f{t) < h{t) for t & I, and 

(A2) the Gegenbauer coefficients satisfy fi >0 for i > t + 1. 

If C G is a spherical r-design of \C\ = N points such that {x,y) G I for distinct 

points x,y & C, then 

(23) E{n,C;h)>N{foN - f{l)). 

In particular, if I = [£{n, N,T),u{n, N,t)], then 

(24) C{n,N,T-,h)>N{foN-f{l)). 
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Proof. Using ([5]) and the conditions of the theorem we consecutively have 


Nf{l)+E{n,C;h) = 


> 


> 


Nf{l)+ 

x,y£C,x^y 

|C|/(1)+ /((!,!/)) 

x,yeC,xjty 

deg(/) f n / \ 2 

iciVo+ E vE 

i=i '■ j=i \xec / 

N^fo, 


which implies that E{n,C]h) > N{fQN — /(!)). If / = [^(n, A^, r), rt(n, r)], the asser¬ 
tion ()24p follows immediately from the definitions (jl6p and (|17p . □ 


Our choice of polynomials for Theorem 13.11 follows from ideas of Levenshtein [331 [3l] 
and the connections cnD and m- We start with fixed dimension n, cardinality N and 
strength r under the assumption that N G {D{n,T), D{n,T + 1)]. Now the equation 
N = LT-{n,s) determines all necessary parameters as explained in subsections 2.3 and 

2.4. 

Definition 3.2. We denote by the set of polynomials satisfying conditions (Al) 

and (A2) of Theorem 13.11 For convenience we shall write An^T-,h '■= 

Next we need Hermite interpolation as follows. If h € 1,1]), we define the 

Hermite interpolant E[t) as follows: 

(i) for odd t = 2k — 1 the polynomial F{t) of degree at most 2 A: — 1 by 

F{ai) = h{ai), F'{ui) = h'{ai), i = 0,1,..., A: - 1; 

(ii) for even t = 2k the polynomial F{t) of degree at most 2k by 

F(/3o) = h(/3o), F(A) = /i(A), T'(A) = h'(A), i = 1,..., A:. 

These conditions define, as in m (see also 121 [Ml m]), a Hermite interpolation problem 
that requires the graph of F{t) to intersect and be tangent to the graph of the potential 
function h{t) at all points ai and all fdi except for /3o = — 1 where only intersection is 
required. 


Lemma 3.3. Suppose h is absolutely monotone on [—1,1] and that F satisfies (i) or (ii). 
Then F{t) < h{t) for all t G [—1,1]. 
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Proof. The proof follows from the well-known error formula for Hermite interpolation 
(see m), namely 


(25) 


h{t) - F{t) 


^(T-+ijP - aof Ok-if, T = 2k-1, 

^ ''(r+i)f^ ~ - I3if • • • (t - /3fc)^ T = 2k, 


for some f, = i{t) G (—1,1) and the fact that > 0 for t G [—1,1], □ 


In |16l Theorem 3.1], the following result was established by the authors, which pro¬ 
vides a lower bound on the /i-energy of general codes. 


Theorem 3.4. Let n,N and r be positive integers with N G {D{n,T), D{n,T + 1)] and 
suppose h is absolutely monotone on [—1, Ij. Then, for any positive integer t, we have 


EiLo T = 2k-1, 


(26) 22(n, N, r; h) > inf E{n,C] h) > < 

c 


,^^E*=o7i^(A), T = 2k, 
where the infimum is over spherical codes C C with ICI = N. 


Utilizing Theorems 13.11 and 13.41 we provide a sufficient condition for the optimality of 
the bounds in ()26p for a class of spherical designs. 


Theorem 3.5. Suppose n, t, and I are as in Theorem \3.1{ N G {D{n,T),D{n,T + 1)], 
h is absolutely monotone on [—1,1] and that ai ^ I for i = 0,... ,k — 1 in the case that 
T = 2k — 1 and that Pi ^ I for i = 0, ... ,k in the case that r = 2k. If C C. is a 
spherical r-design with ICI = N and inner products {x,y) G I for x y £ C, then the 
linear programming lower bounds in ()26[) cannot be improved by utilizing polynomials of 
degree at most r satisfying (Al); i.e., for any such polynomial f we have 


(27) 


N{foN - f{l)) < N{FoN - F{1)) 


Ei=o T = 2k-1, 

^N‘^T,i=olih{Pi), T = 2k, 


where F{t) is the Hermite interpolating polynomial from (i) and (ii). 


Proof. We shall consider only the case r = 2A; — 1 since the t = 2k case is analogous. 
Notice that (i) and (fT^ allow us to rewrite (l26]l as 

k-l 

E{n, C; h)>N^Yl - F(l)). 

i=0 
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Lemma 13.31 implies that F{t) < h(t) for every t G [—1,1]; in particular F{t) satisfies the 
condition (Al) of Theorem 13.11 The condition (A2) is trivially satisfied and therefore 
F G An^2k-l,r,h- 

Furthermore, for any polynomial f{t) G An^ 2 k-i,r,h of degree at most 2k — 1, we have 
from the quadrature formula (fT^ for f(t) and the fact that {«*} C I 

k—1 k—\ 

A^(FoA^ - T(l)) = -/(!)), 

2=0 2=0 

which proves (12711 and the theorem. □ 

In Theorem 14. ll we show that the bound for C{n, N, r; h) in (|26|) can be improved over 
the whole range D(n, t) < N < D[n, r + 1) in the case of even r. 

As mentioned above, the specific properties of the spherical designs, namely the ex¬ 
istence of nontrivial upper bounds on u{n, N,t), allows further application of the linear 
programming techniques. We are able to derive upper bounds on Ll{n, N,t] h) thus 
setting a strip for the energies of the spherical designs under consideration. 

Theorem 3.6. Let n, N, r he positive integers with N > D(n,T) and let h : [—1,1] —S' 
[0, +oo]. Suppose I is a subset of [—1,1) and g{t) = ® polynomial 

sueh that 

(Bl) g{t) > h{t) for t ^ I, and 

(B2) the Gegenbauer coefficients satisfy gi < 0 for i > t + 1. 

If C G is a spherieal r-design of jCj = N points sueh that {x,y) G I for distinct 

points x,y € C, then 

(28) E{n,C-,h)<N{goN-g{l)). 

In particular, if I = [i{n, N,T),u{n, N,t)], then 

(29) U{n,N,T;h)<N{goN-g{l)). 

Proof. Let C C be an arbitrary spherical r-design of jCj = N points. Using ([5]) 

and the conditions of the theorem we consecutively have 

Ng{l) + E{n,C;h) = Ng{l)+ ^ H{x,y)) 

x,y£C,x^y 

< Ng{l)+ ^ 9 {{x,y)) 

x,y£C,x^y 



< N^9o, 
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which implies that E(n,C‘,h) < N{gQN — g{l)). Since the design C was arbitrary, we 
conclude that Z^(n, iV, r;/i) < D 


We utilize Theorem 13.61 to determine an upper bound that in conjunction with the 
lower bound (j27l) in Theorem 13.51 defines a strip where the energy of designs lives. We 
shall formulate the theorem for the odd case t = 2k — 1, but a similar assertion holds for 
the even case t = 2k. 

Theorem 3.7. Let n, N, t = 2k — 1 be positive integers with N G (T)(n, r), D{n, r + 1)] 
and let h : [—1,1] ^ [0,+oo] be absolutely monotone. For a^-i < u < 1 (see (fTT)l ) and 
every j G {0, 1.. . , A: — 1}, let G{t) = Gj^u{t) to be the Hermite interpolant of h of degree 
2k — 1 that satisfies 

G{ai) = h{ai), G'{ai) = h!{ai), i G {0,1,... ,/c -1}\{j}, G(-l) = /i(-l), G{u) = h{u). 
Then, for any spherical r-design G with |G| = N and u{G) = max (x,y) < u, 

x,yGC,x^y 

k-1 

(30) E{n, C; h) < N{GoN - G(l)) = + ^^Pj • 

i=0 

In particular, taking u := u{n,N,T), we obtain 

(31) U{n,N,T‘,h) — C{n,N,T]h) < min pj[Gj^uio(j) — h{aj)]. 

0<j<k—l ’ 


Proof. It follows from the Hermite error formula 

h{t) - G{t) = ^^2kp {t + 1)(A - u) JJ(t - a^)^ 

that G{t) > h{t) for t G [—l,u]. We next apply the quadrature formula (fT3]l to the 
polynomial G{t) and utilize Theorem 13.61 to conclude (1301) . The estimate (l3TI) follows 
from Theorem 13.41 and ([30l) . □ 


4. Improving the linear programming lower bounds 

4.1. Even strength. We show that the bounds from Theorem l3.4l can be improved when 
some additional information about the distribution of the inner products of the designs 
under consideration is available. This is exactly the case for designs of even strength 
2k and cardinality N in the interval {D{n,2k),D{n,2k + 1)). In fact, N = D{n,2k) is 
possible only for k = 1 (the regular simplex) and k = 2 (see [H [5] ). 

It follows from Lemmas l2.3l and [Tn that —1 < i{n, N, 2k) for every N G {D{n, 2k),D{n, 2k+ 
1)). We next use this fact to improve the bound (|26p for the case of even strength r = 2k, 
where N G {D{n, 2k), D{n, 2k + 1)). 
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Theorem 4.1. Suppose N S {D{n,2k), D(n,2k + 1)) and h is absolutely monotone on 
[—1,1]. Let G{t) he the Hermite interpolant of h{t) of degree at most 2k such that 

G{i) = h{e), G(A) = h{^i), G'iPi) = h'iPi), i = l,...,k, 

where i := £(re, A^, 2k'). Then 

k k 

C{n, N, 2k- h) > N{GoN - G(l)) = iV 

i=0 i=0 


Proof. Using Theorem l3.1l with f = G and I = [1, 1], we obtain C{n, N, 2k-, h) > N{GqN— 
G(l)). Since the degree of G{t) is at most 2k we can apply the quadrature formula (fT^ 
and so 


GoiV-G(l) = 


> 


iV^7iG(A) = N ( 7oG(-1) + ^7,/i(/3,) 

i=0 V i=l / 

N (^7o(G(-1) - h{-l)) + ^7,h(/30^ 


k 


i=0 


(the inequality G(—1) > h(—1) follows from the interpolation). 


□ 


We remark that the choice of the polynomial G{t) in Theorem l4.1l is not usually optimal 
for maximizing the lower bound in Theorem 14.11 In a forthcoming paper, we shall develop 
methods for choosing optimal interpolation points. 

Next we show that the inequality for (.{n, N, 2) from Lemma I2.2l fal can be used for 
obtaining a lower bound for C{n,N,2-,h) that is better than (I26p in the whole range 
n + 1 = D{n, 2) < N < 2n = D{n, 3). 


Theorem 4.2. Let n > 2 and N ^ [n + l,2n]. If h is absolutely monotone on [—1,1], 
then 


(32) C{n,N,2-,h)> 
In particular, if C, := N/n, then 

(33) £(IV/C,iV,2; 


N[h{t))N{N - n - 1) + nh{l - N/n)] 
N — n 


h) > h{0)N^ + 


N[h{l - C) - ChjO)] 

C-1 
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Proof. Let n < i{n,N,2). We apply Theorem 13.11 with I = [k, 1] and a second degree 
polynomial f{t) such that /(«;) = /i(k), /(a) = h{a) and f'{a) = h'{a) for some a G (k, 1). 
The optimal choice of a so as to maximize f^N — /(I) is 


(34) 


n(l — k) — N 
n(l — k) + kNu 


By Lemma [T^ ai. £{n, N,2) > 1 — N/n. Choosing k = 1 — N/n gives oq 
the bound (p2]l follows. 


0 and, hence, 
□ 


Lower bounds for the energy of 4-designs by Theorem 13.11 can be obtained by interpo¬ 
lation with polynomials of degree four: 

/(k) = h(K), /(a) = h{a), f'{a) = h'{a), f{b) = h{b), f'{b) = h'{b), 

where k is the expression on the right-hand side of i\22\) and the touching points a and b 
are chosen to maximize JqN — /(I). 

4.2. Odd strength. In contrast to the even strength case, non-trivial lower bounds 
on i(n,N,2k — 1) seem impossible for N G [D(n,2k — 1), D(n,2k)] without further 
constraints. Instead we show that higher degree polynomials can lead to improving the 
lower bound (j26l) for t = 2k — 1. 

Let n, N £ [D{n,2k — 1), D{n,2k)] and t = 2k — 1 he fixed and j be a positive 
integer. Taking the necessary parameters from N = L 2 fe-i(n, s) we consider the following 
functions in n and s G l 2 k-i 

.. k—1 

(35) Qjin,s) := — + '^ p^pj'"\ai). 

1=0 

Note that Qj{n,s) is the same as the quadrature rule expression on the right-hand side 
of (fl^ applied to / = since P^^\l) = 1. 

The functions Qj{n,s) were firstly introduced and investigated in M- The applica¬ 
tions in [T4] target the upper bounds on A{n, s) and, in particular, the possibilities for 
improving the Levenshtein bounds. We are going to see that the functions Qj{n,s) are 
useful for our purposes as well. The next theorem shows that (the signs of) the functions 
Qj{n,s) give necessary and sufficient conditions for existence of improving polynomials 
of higher degrees. 

Theorem 4.3. Assume that h is strictly absolutely monotone. Then the bound i26\) can 
be improved by a polynomial from Anpk-i-,h of degree at least 2k if and only if one has 
Qj{n, s) < 0 for some j > 2k. 

Moreover, if Qj{n, s) < 0 for some j > 2k, then [WH can be improved by a polynomial 
from An^ 2 k-i-,h of degree exactly j. 
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Proof. The necessity mirrors the spherical codes’ analog m Theorem 4.1] (see also |16l 
Theorem 2.6]). We include the simplified proof of the sufficiency in the context of the 
spherical designs for completeness. 

Let us assume Qj{n,s) < 0 for some j > 2k. We prove that the bound (|26]) can be 
improved by using the polynomial f{t) = ePj^\t) + g{t) for suitable e > 0, where the 
polynomial g{t) is such that deg(g') = 2k — 1 and 

(36) g[ai) = h{ai) - £P'f'\ai), g'{ai) = h'{ai) - eiyP'f'^)'{ui), i = 0,1,..., A; - 1 

(i.e. g{t) is Hermite interpolant of h{t) — £P^'^\t) in the points oq, ai,..., afc-i. We 
denote g{t) = giP^'^\t). Note that /o = go and /(I) = ^(1) + e. 

We first prove that f{t) = £P^^\t) + g{t) G for some e > 0. For condition 

(A2) we need to see only that fj = £> 0. For condition (Al), let us choose e > 0 such 

that (ji — ^ (t) >0 for every £ > 0 and for every t G [—1,1]. It is clear that 

such e exists because h is strictly absolutely monotone and this leaves finitely many £ 
to generate inequalities for e. Moreover, since the function h{t) := h{t) — £Pj^\t) is 
absolutely monotone by the choice of e, we could infer as in Lemma 13.31 that g satisfies 
g{t) < h{t) for every t G [—1,1) which implies that f{t) < h{t) for every t G [—1,1] and 
hence f{t) G 

It remains to see that the bound given by f{t) is better than (|26jl . This follows by 
combining equalities from (136]) and applying (fT^ and (l3^ to obtain 

fc-i 

N{Nfo - /(!)) = N^Y. 

i=0 

Since Qj{n, s) < 0, we have N{Nfo — /(!)) > N'^ Yi=o i.e. the polynomial f{t) 

indeed gives a better bound. □ 

As mentioned above, the test functions Qj{n, s) were initially defined in [T3] as related 
to the Levenshtein bounds C{n, s) on maximal spherical codes. Theorem 14.31 shows 
that Qj{n,s) prove to be very useful tool in the context of potential energy as well. In 
particular, the signs of the test functions Q 2 k+ 3 {n, s), where r = 2A; — 1, were investigated 
in detail in [T3] (see also m)- Denote 

, - 4A; + 5 + ^/W^^W^^6k^T2AkT^ 

ko--= - 4 - 

for short (fco is well defined for k >9). 

Theorem 4.4. [TO] IFe have Q 2 k+ 3 {n,s) < 0 for every s G 
n > 3 and k >9 which satisfy 3 < n < ko. 


and for every 
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Corollary 4.5. The bound ![26\) can be improved by using polynomials of degree 2k + 3 
for every s € /o’" every n > 3 and k > 9 that satisfy 3 < n < ko. 

Proof. This follows from Theorems 14.31 and 14.41 □ 


Remark 4.6. We note that test functions for the even case t = 2k can be defined and 
investigated as well (see M)- However, we conjecture that the linear programming in 
[£(n,N,T),u{n,N,T)] by Theorem 13.11 is always better than the bounds which would 
come from higher degree polynomials. 


5. Improved upper bounds for 2, 3, and 4-designs 

In this section, we use bounds from Lemmas o and 12.21 to specify an interval I in 
Theorem l3.6l to obtain upper bounds for the energy of 2-, 3- and 4-designs. Our numerical 
experiments suggest that the use of even degree polynomials is not effective and so we 
turn our attention to polynomials of degrees 1 and 3. 

5.1. Upper bounds for 2-designs. First we apply Theorem 13.61 for g a linear polyno¬ 
mial. 


Theorem 5.1. Let h be a convex non-negative function on [—1,1] and let u and i denote 
the upper and lower bounds in dii) and (EH), respectively. For N £ [n -|- 1, 2n), if i < u 
then 


(37) 
ifi = 


U{n, N, 2; h) < 


N[{N - l){uh{i) - £h{u)) + h{e) - h{u)] 
u — i 


u then U{n, N, 2; h) = Nh{—1/{N — 1)). 


Proof. With I = [f, n], the linear polynomial passing through the points {i,h{i)) and 
{u,h{u)) satisfies the conditions of Theorem 13.61 and gives the desired bound. If n = 
i, then we must also have u = i = —1/{N — 1) which implies that U{n, N,2;h) = 
Nh{-1/{N -1)). □ 


Remark 5.2. Combining (1321) and (1371) gives a strip for the /i-energy of any spherical 
2-design of G {n -|- 1, n -|- 2,..., 2n — 1} points when h is absolutely monotone. Note 
that if n and N tend simultaneously to inhnity such that n/N —)• ^ for some C £ (1,2), 
Theorem 14.21 and Theorem EU give an asymptotic strip as N ^ oo: 


( 38 ) 


MI-O^ + o(iv-) < 

/I /-\AT \ / J\J-Z 


< 


(I-ON 

hji-O + hjC-i) (2-C)h(i-C)-CMC-i) , ^ ( ^-2^ 

2(C-1)A^ 


2 
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Example 5.3. Simple algebraic manipulations show that the bounds (1321) and (I37p for 
2-designs coincide when N = n + loi:N = n + 2 for every n and h (i.e. the strip becomes 
a point for these two cardinalities). The case = n -|- 1 leads to the regular simplex on 

The case N = n + 2 is more interesting - Mimura [36] has proved that spherical 2- 
designs with n -I- 2 points on do exist if and only if n is even and Sali m (see also 
|38j i proved that there are no other (up to isometry) such 2-designs. Spherical 2-designs 
of = 2A: points on are known as Mimura spherical designs and consist of two 

orthogonal fc-simplices which we denote by {A:, k}. This design has k{k — 1) distances of 
y^2k/{k — 1) coming from edges within the two /c-simplices and k"^ distances of \/2 which 
are the edges joining the vertices from distinct simplices. The total number of various 
distances is k{2k — 1) = ('^). 

Sali’s nonexistence result follows easily from the coincidence of our bounds. It also 
follows that the 2-designs of n -|- 2 points for even n are unique and optimal - they 
have simultaneously minimum and maximum possible energy. The optimality cannot be 
extended to the larger class of spherical codes - Cohn and Kumar m Proposition 1.4] 
prove that if n-|-l < N < 2n, then there is no A^-point universally optimal spherical codes 
on For the cases, N ^ {n -|- 3, n -|- 4,..., 2n — 1}, there is a difference (increasing 

with N) between the bounds from (|32p and (I37p . 


Example 5.4 (Riesz s-energy of Mimura designs). As mentioned above, Mimura designs 
clearly give the minimum energy over 2-designs with 2k points in The Riesz s- 

energy, that is the energy when h{t) = (2(1 — t))“^/^, is given by 


(39) 


Esiikk}) 


e / k{k-l) fk-iy^‘^\ 

^ ^ P \ k ) j’ 


It is easy to see that Mimura designs are not universally optimal among general codes 
of cardinality N = 2k. Indeed, consider the competing configuration {2, 2k — 2} made of 
two orthogonal simplices, a 2-simplex and a {2k — 2)-simplex (say the North Pole-South 
Pole diameter and a (2A; — 2)-simplex in the equatorial hyperplane). This conhguration 
has 1 distance of length 2, {2k — 2){2k — 3)/2 distances of length ■\/2{2k — 2)j{2k — 3) 
and 2{2k — 2) distances of length y/2. Thus the s-energy of this competing conhguration 
is given by 


(40) 


Es{{2,2k-2}) 


2{2k - 2) / 1 2k-3 f 2k - 3 

2^/2 2{2k - 2)2^/2 ^ 4 \2k-2) 


thus showing that for large enough s the Mimura design {k, k} will have strictly larger 
energy than the competing conhguration {2,2k — 2}. However, for k = 3 it was shown 
in |22| that the Mimura conhguration minimizes, in particular, the logarithmic energy. 
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5.2. Upper bounds for 3 and 4-designs. For our estimates we shall apply Theo- 
rem l3.6l with g{t) the Hermite interpolating polynomial of degree at most three satisfying 

(41) g{£) = h{e), g{a) = h{a), g'{a) = h'{a), g{u) = h{u), 

where i and u again denote lower and upper bounds, respectively, for the inner products 
of all spherical 3- or 4-designs under consideration. For r = 3 one can take i = —1 and 
u as in [131 Theorem 3.9]. For r = 4 the bounds u and £ are taken from Lemma l2.1f bi 
and l2.2f bL respectively. 


Theorem 5.5. Let h G C'^([—1,1]) with h^^^t) > 0 for t G [—1,1]- For r = 3 and 
N G [2n, and for r = 4 and N G -|- n), we have 

(42) U{n,N,T;h) < N{N-l)h{ao) 

{h{£) — h{ao)) [rtA^(l -|- nog) -|- 2Nao -|- n(l — «)(! — oq)^] 

^ n{u -£){£- ao)^ 

{h{u) - h{ao)) [£N{1 + nug) + 2A^ao + n(l - £){l - uq)^] 

^ n{u-£){u-aof ’ 


where 


(43) 


N{1 -|- u) -|- ra(l — 1){1 — u) 
n(l — £){1 — u) — N{1 + £un) ’ 


and u and £ are chosen as described above. 


Proof. As in Lemmawe see that g{t) > h{t) for every t G [£,u], i.e. the condition 
(Bl) is satisfied. Moreover, (B2) is trivially satisfied and therefore Theorem 13.61 can be 
applied. Standard calculations for optimization of g^N — g{l) via derivatives show that 
the optimal value of a is given by (f43l) . □ 

The asymptotic form of the bound (|42[) for r = 4 is easily determined when n and 
N G [Zl(n, 4), Zl(n, 5)) = pU+T ^ ^^2 _j_ tend simultaneously to infinity as described in 
the following corollary. 

Corollary 5.6. Ifn and N tend to infinity in relation N = n^A-l-o(l) as N ^ oo, where 
Ag[1/2,1) is a constant, then 

(44) U{n,N,4:-,h) <h{0)N^-h{0)N+ ciVn+ C 2 + o{l), {N ^ oo), 

where 

_ V A[(2Va - l)h{l - 2VA) + (1 - 2VX)h{2Vx - 1)] 

2(2VA - 1)3 

(1 - VX)h{l - 2^/X) + Va/i(2Va - 1) - h{0) 

(2VA-1)3 ■ 


and 
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Proof. The assertion follows from Theorem 15.51 and the asymptotic formulas below: 
u{n,N,4) = 2VA-1 + o(1) (from Lemma l2.11 bll. 

£{n,N,4:) = l-2VX + o{l) (from Lemma l2.2f bll. 
ao = o(l) (from (Ii3])l. 


□ 


6. Some asymptotic lower bounds 


We consider the bounds (1261) in the asymptotic process where the strength r is fixed, 
and the dimension n and the cardinality N tend simultaneously to infinity in the relation 


(45) 


lim 


N 


+ 7 


n,N->-oo ^ {k — 1)! 

(here 7 > 0 is a constant and the term comes from the Delsarte-Goethals-Seidel 

bound). 

Theorem 6.1. Let h he absolutely monotone on [—1,1] and r be fixed. If n and N tend 
to infinity as in (SSI), then 

(46) C{n,N,T-h)>h{{))N‘^ + o{N'^). 

Proof. Let t = 2k — 1. From ()45l) and |12] . the asymptotic behavior of the parameters 
from (f26]l is as follows: Oj = o(l), for i = 1, 2,... , A: — 1, ao = — + o(l)) 

PqN = (1 + 'y{k — 1)!)^^“^ + 0(1) as —)• 00. Now the lower bound of (126|) is easily 

calculated: 

k-l 

C{n, N,2k — l;h) > N'^ E pMai) 

i=0 

= ^Po/i(ao) + h{^) ^ + o{N‘^) 

= (^po{h{ao) - hm + h{0) (^1 - ^ + o(iV2) 

= h{0)N^+ o{N'^). 

Similarly, in the even case r = 2k, we obtain using the appropriate asymptotic relations 
from [12], 

C{n,N,2k-h) > + 


= hfi))N^ + o{N^) 


□ 
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Note that for r = 4 the main terms in the lower bound (j46p and the upper bound 
([441) coincide. For odd r, (f46l) also yields good energy estimates as shown in the next 
example. 


Example 6.2. There is standard construction (see m Chapter 5]) mapping binary 
codes from the Hamming space H{n, 2) to the sphere - the coordinates 0 and 1 are 
replaced by ±l/^/n, respectively. Denote by x and C the images of vector x € H{n,2) 
and code C C H{n, 2). Then the inner product {x,y) on and the Hamming distance 
dH{x,y) in H{n,2) are connected by 


{x,y) = 1 - 


2dH{x,y) 

n 


Levenshtein [33l pages 67-68] shows that spherical codes which are obtained by this 
construction from the Kerdock codes [28] are asymptotically optimal with respect to their 
cardinality. So it is natural to check them for energy optimality. 

The Kerdock codes Ki C H{2‘^^,2) are nonlinear. They exist in dimensions n = 2^^ 
and their cardinality is = 2^^. The Hamming distance (weight) distribution does 

not depend on the point and is as follows: 

Ai = 0 for i 0,ii = 2^^-^- 2 ^-\i 2 = 2 ^^-\i 3 = 2^^-^+ 2^-\n = 2^^, 

^0 = 1; 

A, = 2^^(2^^-^-l), 

= 22^+1-2, 

= 22 ^( 22 ^-! - 1 ), 

An = 1. 

The weights 0, H, ^2, is, u correspond to the inner products 1, 0, —— 1, respectively. 

The spherical code K£ C §2^^-i has energy 

E(n,K,;/i) = iv(^(22^+i-2) h(0) + 2^^(2^^-^ - 1) (^h + h +M-1)) • 

When n tends to infinity we obtain 

_ j \ t 2 /I 1 \ 

E(n, Ke; h) = — h(^) + h( -+ 0(N) = h(0)N^ + 0(N^/^) = h{0)n^ + 0{n^), 

2 \ / 

where we assumed that h is differentiable at 0. 
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